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Abstract 

A new approach to probability theory based on quantum mechanical and Lie algebraic 
ideas is proposed and developed. The underlying fact is the observation that the co¬ 
herent states of the Heisenberg-Weyl, su{2), su{r-\- 1), sri(l, 1) and su{r, 1) algebras in 
certain symmetric (bosonic) representations give the “probability amplitudes” (or the 
“square roots”) of the well-known Poisson, binomial, multinomial, negative binomial 
and negative multinomial distributions in probability theory. New probability distri¬ 
butions are derived based on coherent states of the classical algebras B^, Cr and Dr 
in symmetric representations. These new probability distributions are simple general¬ 
isation of the multinomial distributions with some added new features reflecting the 
quantum and Lie algebraic construction. As byproducts, simple proofs and interpre¬ 
tation of addition theorems of Hermite polynomials are obtained from the ‘coordinate’ 
representation of the (negative) multinomial states. In other words, these addition 
theorems are higher rank counterparts of the well-known generating function of Her¬ 
mite polynomials, which is essentially the ‘coordinate’ representation of the ordinary 
(Heisenberg-Weyl) coherent state. 


^On leave of absence from Institute of Theoretical Physics, Northeast Normal University, Changchun 
130024, P.R.China. 



1 Introduction 


Quantum theory is one of the greatest achievements in twentieth century physics. It has 
changed the fundamental structure of physics, material science and also influenced various 
disciplines, in particular biological (genetic) science and philosophy. Quantum theory dic¬ 
tates that at the microscopic level nature is not governed by causal laws typically exemplified 
by the Newtonian equation of motion but by probabilistic laws. The fundamental ingredient 
of quantum theory is, however, not the probability itself but the probability amplitude which 
obeys a certain equation of motion and the square of which gives appropriate probabilities. 

In the present paper we report on an attempt to apply quantum theory ideas to probabil¬ 
ity theory itself. This, we believe, will provide new perspectives on probability theory and 
hopefully will enrich the long-established and rather mature science. The hrst step would 
be to associate certain “probability amplitudes” to some typical probability distributions of 
classical probability theory. In a broader perspective, this problem belongs to the paradigm 
of “square roots”. The Dirac equation is obtained as a “square root” of the Klein-Gordon 
equation. The creation and annihilation operators can be considered as “square roots” of 
the harmonic oscillator hamiltonian. Of course such a “square root” can never be unique. 
It depends on the formulation. It turns out that the ^coherent states^ B I, I, i in quan¬ 
tum optics and the so-called ‘generalised coherent states’ HB associated with various Lie 
algebras could be identihed as certain “probability amplitudes”. For example, the coher¬ 
ent states associated with the Heisenberg-Weyl algebra, su{2) |], su{r -|- 1) [T^ and 

su(l, 1) B H, |TT|, 1^, 1^ su{r, 1) [T^ algebras in totally symmetric (bosonic) representations 
could well be interpreted as “probability amplitudes” for the Poisson, binomial, multinomial 
and negative binomial, negative multinomial distributions in probability theory, respectively 
131, IT3. This also means, in turn, that these typical discrete probability distributions are 


characterised in terms of Lie algebras (groups) and their representations. The relationship 
between the Poisson distribution and the ordinary coherent states is well-known and that 
of the binomial distribution and the 5^(2) coherent states is also known, but to a lesser 
degree. The characterisation of the negative binomial (multinomial) distributions by Lie- 
algebra representations, which is believed to be new, has been reported in our previous work 


The second step is to extract useful information (predictions) from the characterisation 
“probability amplitudes = coherent states”. One would naturally ask ‘what would be the 
probability distributions associated with the other Lie algebras and/or other representa- 
^In this paper we call them simply coherent states. 
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tions?’ In the present paper we mainly address the problems in this step. We choose the 
classical Lie algebras, Bj., Cr and in Cartan notation (or so(2r + 1), sp{2r) and so(2r) 
algebra, respectively) and construct the coherent states in the totally symmetric (bosonic) 
representations. This gives rise to new probability distributions, to be denoted as Bj. multi¬ 
nomial distributions, etc. One reason for choosing the symmetric representations is that 
they are supposed to give closest analogs of the classical probability distributions, like the 
multinomial distribution. Another reason is the relative ease of the calculation and presen¬ 
tation. 

The third step would be to discuss the time evolution (stochastic process) based not on 


the probability itself but on the “probability amplitude” in the spirit of quantum theory 
This would be the subject of our future publication. 


This paper is organised as follows. In section two we explain the basic idea of intro¬ 
ducing the “probability amplitude” by taking the simplest and well-known example of the 
Poisson distribution and derive the ordinary coherent state. This section is meant for wider 
readership. In section three we discuss the “probability amplitudes” for the binomial and 
multinomial distributions, the coherent states of Ai (su(2)) and Ar {su{r -|- 1)) algebras in a 
slightly different way from our previous work [Q . The representation theory aspects of these 
algebras are emphasised in order to facilitate the transition to the other algebras treated in 
later sections. As new material in this section we discuss the x (coordinate) representation of 
these coherent states. Based on new expressions of the Ai and Ar coherent states, which have 
straightforward interpretations of “probability amplitudes” for the binomial and multinomial 
distributions, we obtain a simple (quantum theoretical) proof and interpretation of addition 
theorems of the Hermite polynomials describing the number states of harmonic oscillators. 
This is analogous to the well-known fact that the coordinate representation of the coherent 
state of the Heisenberg-Weyl group gives the generating function of Hermite polynomials. 
In sections four, five and six, we derive new probability distributions associated with the to¬ 
tally symmetric (bosonic) representations of the Cr, Br and Dr algebras, respectively. These 
are the first and simplest results of the second step of the “quantum theory of probability” 
mentioned above. Since the Dynkin diagram of Cr is obtained from that of A 2 r-i by folding, 
the Cr coherent states resemble closely those of the A 2 r-i algebra. However, the obtained 
probability distributions, to be denoted as the Cr multinomial distributions, have markedly 
different features from the ordinary multinomial distributions, reflecting the different weight 
space structures of the Cr and A 2 r-i algebras. The probability distributions associated with 
the symmetric representations of Br and Dr algebras have also new and interesting features. 
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Since Br Dynkin diagram is obtained from that of -D^+i by folding, these probability dis¬ 
tributions are somewhat related. Section seven is devoted to a summary of results. In the 
Appendix we give a simple proof and interpretation of another type of addition theorems of 
Hermite polynomials based on the x representation of sm( 1, 1) and su{r, 1) coherent states. 
The formula is known as generalised Mehler formula but is not found in the standard mathe¬ 
matics reference texts. This time the summation includes inhnite number of terms reflecting 
the inhnite dimensionality of the irreducible unitary representations of these non-compact 
algebras. 

2 “Quantum Theory of Probability”: An Example 

Let us begin with the naive idea of associating “probability amplitude” to a probability 
distribution. In other words, we explain how to give some meaning to a “square root” 
of a probability distribution by taking the simplest example of the Poisson distribution. 
Throughout this paper we consider only discrete probability distributions P parametrised 
by a set of integers. A probability distribution parametrised by one non-negative integer n is 
completely specihed by a set of non-negative numbers satisfying the conditions of unit total 
probability: 

OO 

Pn>0, E^n = l. (2.1) 

n=0 

For a quantum theory let us introduce a Hilbert space Ti with an orthonormal basis |n), 
n = 0,1,2,..., 

{tTi\ti) ^mni ( 2 - 2 ) 

satisfying the completeness relation 

OO 

( 2 - 3 ) 

n=0 

in which I on the left hand side is the identity operator. Our objective is to hud a nor¬ 
malised state 1'^) in Ti such that its transition amplitudes (n|'^) give rise to the probability 
distribution: 

|(n|^)|2 = P„, n = 0,1,2,.... (2.4) 

Then by using the completeness relation one obtains 

OO OO _ 

1^) = I] \n){n\il)) = Y, (^.5) 

n=0 n=0 

in which the phase 5n is arbitrary. Thus far the Hilbert space is unspecihed. 
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Let us choose as Ti the Hilbert space of one of the simplest quantum systems, the har¬ 
monic oscillator. It is described by the annihilation and creation operators a and satisfying 
the commutation relation 

[a, a^] = 1. (2.6) 

(Throughout this paper Planck’s constant h is set to unity.) Then the orthonormal basis is 

simply given by 

|n) = ^4|0), n = 0,1,2,..., (2.7) 

\/nl 

in which |0) is the vacuum state characterised by the condition 

a|0) = 0. (2.8) 

The well-known Poisson distribution describing random processes occurring in a time 
(space) sequence is 

2 

F„(a)=e-“—, n = 0,1,2,.... (2.9) 

nl 

For example, the number of radio-active decay particles emitted from a sample in a hxed time 
(t) is known to obey this distribution, oc t. Then the quantum state |'0(a)) (“probability 
amplitude”) corresponding to the Poisson distribution (|2.9|) is easily obtained (we set = 0): 

oo 

W.(a))=e-”’'''^i:4=|«). (2.10) 

.1.0 V 

If we substitute the dehnition of the number state in terms of the creation operator, we 
obtain a closed form 

\^(a)) = = e“(“^-“)|0), (2.11) 

and the last formula is obtained by using the Baker-Campbell-Hausdorff (B-C-H) formula 

gA+S ^ gAgBg-i[A,i?] 

for the case [A, B] commutes with A and B. This state was hrst introduced by Schrodinger 
1^ and discussed by many authors ^ under the name ‘coherent state’ which was coined 
by Glauber in quantum optics. The coherent state has many other characterisations. 

1. It is an eigenstate of the annihilation operator: 

01-0(0:)) = o| 0 (o)). 

2. It is a minimum uncertainty state: 

{Ax‘^){Ap^) = 1/4. 
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in which x = (a^' + a)/\/2, p = i(a^ — a)/\/2 are the corresponding coordinate and 
momentum of the oscillator. Heisenberg’s uncertainty principle dictates that 

(Aa:^)(Ap^) > 1/4, 


for arbitrary states. 

3. It is obtained by applying a unitary operator (known as the displacement operator) 

ga(at-a) 


to the vacuum state. Such unitary operators form a (unitary) representation of the 
Heisenberg-Weyl group. 


The last characterisation is generalised by many authors and the concept of the coherent 
states associated with various Lie algebras (groups) is now well established. Thus starting 
from a rather naive idea of introducing “probability amplitude” for the Poisson distribution 
we have arrived at the concept of the coherent states, a rather solid subject in quantum 
theory and the representation theory of Lie algebras (groups). As we have shown in previ¬ 
ous publications O, M, the relationship between coherent states and certain probability 


amplitudes is neither coincidental nor superficial but essential. As we will briefly review in 
the next section, the “probability amplitudes” for the well-known binomial and multinomial 
distributions are the coherent states of su{2) and su{r + 1) algebras in the totally symmetric 
(bosonic) representations. The same assertion holds for the negative binomial and nega¬ 
tive multinomial distributions and the corresponding algebras are su{l, 1) and su{r, 1), the 
non-compact counterparts of su{2) and su{r + 1). 


3 Coherent States of Ar algebra 

3.1 Binomial States 

Let us continue along the line of argument of introducing “probability amplitudes” for clas¬ 
sical probability distributions. Here we consider the binomial distribution: 

B{no,nAv'^M) = no + ni = M, peK, (3.1) 

which describes probability distribution of M Bernoulli trials of success (probability p"^) and 
failure (probability 1 — p"^). Here rii is the number of successes and no failures. As a Hilbert 
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space let us choose the Fock space generated by two independent bosonic oscillators: 

[ttj, O/j] Sjk, [oj, Ofc] [oj, O/j] 0, j, k 0,1, 

Km) = a,|0) = 0, j = 0,l, (3.2) 

and restrict the total number to M (integer) 

no + rii = M. (3.3) 

Let us denote by \r]] M) the “square root” of the binomial distribution within this finite 
(M + 1) dimensional Hilbert space. Following the same steps as in the previous section, we 
arrive at a simple expression: 

\mM) = 

nQ-\-n\=M 



which shows clearly that the “transition amplitude” for each possible result (no, Uil?]; M) is 
actually obtained by the binomial expansion. 

The next step is to identify |? 7 ; M) as a coherent state. Let us recall the realisation of 
su{2) algebra in terms of two bosonic oscillators: 

J+ = aoUi, J- = a|ao, Jo = -(^o^o ~ o|ai), 

[J+, J-] = 2Jo, [Jo,J±] = ±J±. (3.5) 

Obviously the restricted two boson Fock space provides the irreducible (spin M/2) represen¬ 
tation of su{2) corresponding to the Young diagram 

I I I I ■ ■ ■ m M boxes. 

Its normalised highest weight state is 

= -^+1^.0) = 0, J„|M,0> = ^|M,0). (3.6) 

Similarly to the coherent states of the Heisenberg-Weyl group in the previous section, su{2) 
coherent states have the form 

H|^o), UeSU{2). (3.7) 
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These coherent states have “minimal uncertainty” if the ‘base’ state |V’o) corresponds to 
a dominant weight, i.e., to the highest weight state or its trajectory by the Weyl group 
|T^ . Thus without loss of generality we choose IV’o) = 1^)0). Since J+ annihilates the 
highest weight state and Jq does not change it, the non-trivial action is by J_ only. So the 
un-normalised su{2) coherent state is given by 

1 ^ IV 1 


e«'^-|M,0) = 


=ef“l«(4)"|0> = -^(4 + ^al)"|0), ? € C, 


(3.8) 


/M! VM! 

Here use is made of the fact that the oscillator algebra [oq, Oq] = 1 is realised by Oq = 8/dal 
and aj). At the last equality, the formal Taylor’s theorem 

e^^f{x) = f{x + a) (3.9) 


is used. It is easy to get the normalised coherent state 

(\/l- \v?al + ria\j |0), r] = ^/+ e C, (3.10) 

which has the same form as the binomial state derived above. (In order to get complex r] we 
only have to choose the phase of M) appropriately.) Thus we have shown that 

the “probability amplitude” of the binomial distribution is the su{2) coherent state. 


3.2 Multinomial States 


In this subsection we discuss the relationship between the multinomial distributions and the 
Ar coherent states B, which has been demonstrated in some detail in our previous paper 


14|. Here we give a simpler and clearer proof of the correspondence with more emphasis on 


the Lie algebraic structures (i.e., roots and weights) which would be useful for comparison 
with the results of the other algebras discussed in later sections. 


The multinomial distribution is 

M' 

Mn(?7;M) = ——^—7 no + niH-h = M, (3.11) 

rif)'. ■ ■ ■ Tlr'. 

in which 

n= = 0 < 17^ = 17J H-Hi?? < 1, e R, J=0,...,r. 

(3.12) 

As a Hilbert space let us choose the Fock space generated by r + 1 independent bosonic 
oscillators 

[aj,a\] = 6jk, aj|0) = 0, j = 0, l,...,r, 

= ^|0), (at)"= (4)"°(4)"i---(4)-y n! = no!ni!---nd, (3.13) 

Vn! 
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n) 














and restrict the total number to be M 


rio + rii + ■ ■ ■ + rir = M. 


(3.14) 


(3.15) 


It has the dimension 

/ M + r\ /M + 

[ M ) = [ r , 

Let us denote by [ 77 ; M) the “square root” of the multinomial distribution within this Hilbert 
space. Then we obtain in a similar way to the binomial state 

HM) = \no, - ■ ■ ,nr){no, - ■ ■ ,nr\'n]M) 

liqH- \-nr=M 

^ \/M\ 


\/no\ ■ ■ - rirl 




M\ 


(hoaSr---(Wrio) 


M\ ^ no!ni! ■ ■ - n,.! 

(?7oaJ) + ?7 i 4 H-h?7r4) |0). 


'M\ 


(3.16) 


Now let us consider algebra and its representations. Its Dynkin diagram is a simple 
line connecting r vertices. The number attached to each vertex corresponds to the name of 
the simple roots given below. 

2 r—1 r 


1 

o 


-o 


The simple roots are most conveniently expressed in terms of r + 1 orthonormal vectors in 

Cj ■ Cfc = 5jk, j,k = 0,l,...,r: 

Oil ~ ^0 — ^1: 0(2 — Cl — C 2 , ■ ■ ■ ) Oir ~ C)—1 — Cj.. (3.17) 

Then any root, positive or negative, can be expressed as 

ej - Cfc, j ^ k, (3.18) 

which is positive ii j < k and negative for j > k. All the roots have the same length. The 
fundamental weight vectors, {\j]j = 1,... ,r}, the dual basis of the simple root system 

2Xj ■ ak/al = 6jk, (3.19) 

can also be expressed by {cj}. For example 
1 

Ai = pY tti + (r - 1 ) 0:2 H-h o,,) = Co - (co + Cl H-h e,.) /(r + 1). (3.20) 


9 



We consider the irreducible representation of Ar with the highest weight 

p = MXi = Mcq — M (eg + Cl + • • • + e,.) /(r + 1), (3.21) 

corresponding to the Young diagram 

rrn ---m m boxes, 

which has the same dimension as the restricted multiboson Fock space introduced 

above. Thus this completely symmetric representation can be realised in terms of r + 1 
bosonic oscillators. The weights and the occupation numbers are related one to one, namely 
the state |no, rii,..., rir) has the weight 

r 

/r = ^ njCj - M (eg + Cl H-h e^) /(r + 1). (3.22) 

i=o 

All the weight spaces are non-degenerate, i.e., one-dimensional. 

If we denote the Ar generators corresponding to the root Cj — Ck by we have 

^U-k) = OjOfe (3.23) 

and 

[^{j-k),X{k-l)] = [o-]o-k,0.iO‘l] = 

[Y(j , Y(^ _j)] = cijCij (3.24) 

Here belongs to the Cartan subalgebra. The quadratic Casimir operator is 

V 

C 2 = — —Xtot{Xtot + r -|- 1), Wot = (3.25) 

which takes the value rM{M -|- r -|- l)/(r -|-1) in the present representation. The state having 
the highest weight (|3.21|) is 

|M,0,...,0) = ^|0), (3.26) 

which is annihilated by the generators 

J, = 1, ■ ■ ■, A (3.27) 

forming an Ar-i subalgebra. The action of the Cartan subalgebra generators 77(0.1) does not 
change the state, either: 


77(o,,)|iV7,0,...,0)=M|M,0,...,0). 
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Thus the coherent states based on the highest weight state ( p.21|) are characterised by 


SU{r + l)/t/(l) X SU{r) = CPh 


(3.28) 


Among the generators belonging to CP*^, only those 


^(y— 0 ) i I 5 • • •! ^ 


(3.29) 


have non-trivial action on the highest weight state ( p.21j) . Thus we hnd, as in the case of 


the binomial state (|3.8|) , that the un-normalised Ar coherent state is expressed as 


e&i«^“])“o(at)^|0) 


'M\ 


M 




+ |0)’ 


| = (6,■■■,e.)eCP^ 


(3.30) 


i=i 


in which use has been made of the Taylor expansion theorem (p.9|) with uq = d/da^. 


The normalised A^. coherent state in the totally symmetric representation is given by 

\n\M) =-j= U)oal + Y.ViA |0). % ={yvlT|ip 6 C, t)o = \/l - |t7P, 

(3.31) 

which has the same form as the multinomial state [ 77 ; M) derived above. As in the binomial 
state case the “transition amplitude” (uq, ..., M) to each number state (or weight state 
(/ii,..., /ir| 77 ; M)) is simply obtained by multinomial expansion. 


3.3 Coordinate Representation and Addition Theorems of Her- 
mite Polynomials I 


In this subsection we consider the ‘coordinate representation’ of the multinomial state ( |3.31| ). 
This representation is useful in quantum optics. It also gives a simple proof and interpretation 


of the following addition theorem of Hermite polynomials (see, for example, [jT^ and pl96 
of m): 


+ ■ ■ ■ + VrXr)l^ril + --- + Vpj 

E (3.32) 

n,+.^nr=M nr\ 
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Here t]o,. .. ,7]^ are arbitrary complex numbers. It should be noted that the left hand side 
contains \/Vo Vr even powers only, since Hermite polynomials have a dehnite parity: 

Let us begin with a single boson oscillator 


[a, a^] = 1. 


The coordinate representation of the number state \n) is 

(xln) = = — - —— =Hn{x) , 

^ ^ IV ; I / ^l/42n/2y^ ^ 

in which Hermite polynomial Hn is given by Rodrigues formula: 

H4x) = {-ire^^D-e-^\ D = 

ax 

It is well-known that the generating function of the Hermite polynomials 

oo 


E ^Hn{x) = e 

n=0 




(3.33) 


(3.34) 


(3.35) 


is essentially the same as the coordinate representation of the coherent state of the Heisenberg- 
Weyl group ( |2.1CI|) : 


{x\^jJ{a)) = e 2 


- ^-Ux-V2a)^ l^llA 


jir' , a G R. 


(3.36) 


The coordinate representation of the multinomial state ( 3.31 ) is simply obtained by expansion 
{Vi, ... ,Vr are in general complex): 


{xo,xi,.. .,Xr\'n-, M) 

-.{X(i,Xi,...,Xr\ (voal^ -h?7r-4)^|0) 


'M\ 

,_ p-ki^O + '-'+^r) 

= VW. .. 

rqH- \-nr=M 


70 7r 


7rh+i)/42^/2 ^ . no! 


HriQ (^0 ) ' ' ' i^Xr ). 


(3.37) 


Next we consider operators A and A dehned by 


A = 


Voo-o 


Vrttr 


ho + ■ ■ ■ + hr 


A = 


Voal H-h Vral 


ho + ■ ■ ■ + hr 


(3.38) 


They are not hermitian conjugate of each other but they satisfy the same relations as those 
of the single oscillator: 

[H,H] = 1, H|0)=0, 


12 














which are essential for deriving Herniite polynomials. Thus we obtain 


{xo,Xi, ... ,x^|77;M) 

(ho + ■ ■ ■ + hr 


{Xo,Xi, . . .,Xr\A^\Q) 




(ho + ■ ■ ■ + hr)^'^^ 6 


H 


M 


{r]oXo H - h r]rXr)/ - ^ r]f^ . (3.39) 


7rh+i)/42^/2 

Comparing (|3.37|) and (|3.39|) we obtain the above mentioned addition theorem ( p.32|) of 
Hermite polynomials, which is nothing but the multinomial expansion of the multinomial 
state. In the Appendix we give a proof and interpretation of another type of addition 
theorems of Hermite polynomials based on negative multinomial states, i.e., the coherent 
states of su{r, 1) algebra in discrete symmetric representations. 


4 Cr Multinomial States 

Let us proceed to the second step in the study of “quantum probability”. In the previ¬ 
ous sections we have shown that some of the typical discrete probability distributions are 
characterised by Lie algebras through coherent states. Now we reverse the logic and try 
to derive new probability distributions starting from Lie algebras and their representations. 
For this we have, in principle, an inhnite choice of Lie algebras and their representations. 
Probably most of such new probability distributions are too exotic to have any practical use 
at the moment. However, the great role played by the Poisson, the binomial, the multino¬ 
mial distributions and their “negative” (non-compact) counterparts makes us expect that 
the probability distributions related with the totally symmetric representations of the other 
classical algebras, Br, Cr and Dr could be useful, though possibly to a lesser degree. Apart 
from the Poisson distribution which has only one parameter, the (negative) multinomial dis¬ 
tribution has many parameters, rj and M, to give suitable description to various statistical 
phenomena. The same property is shared by all the probability distributions derived from 
the totally symmetric representations of Br, Cr and Dr algebras. We propose to call these 
coherent states the Br, Cr and Dr multinomial states and the corresponding probability 
distributions the Br, Cr and Dr multinomial distributions. We start with the Cr case and 
proceed to Dr and Br cases, in the order of increasing complexity. 

4.1 Coherent States 

The Dynkin diagram of C is obtained from that of A 2 r-i by folding. 
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2 


2 


1 

o 


r 



r—1 



1 

a 



r 



2r-2 2r-l 



Its simple roots can be expressed most conveniently in terms of an orthonormal basis of R*”, 

‘ jjk 0,..., r. 

Oil = ei — 62 , 0!2 = 62 — 63 , • • • , CKj—l = — Cj., CKf = 2er. (4-1) 

The positive roots are 

ej-Ck, {j<k), ej + Ck, 2ej. (4.2) 

There are 2r(r — 1) short roots and 2r long roots (±2ej) and the dimensions of Cr algebra 
is 2r^ + r. The fundamental weights are 

Ai = Cl, A 2 = Cl + 62 , ... (4-3) 

We consider the irreducible representation with the highest weight 

p = MAi = Mei. (4.4) 


Its dimensionality is 

/M + 2r - l\ _ /M + 2r - l\ 

1^2r-lj“l^ M ;■ 

It is the same as the dimension of the restricted multiboson (M particle) Fock space of A 2 r-i 
with 2r bosonic oscillators: 


K') «I] = Ikj, bl] = Sjk, j,k = l,...,r 

with the number states 


(4.5) 


\ni,... ,nr-,ni,... ,nr), rii -h + ni H-h = M, (4.6) 

in which Uj {fij) is the number of Qj {bj) quanta. 

Similarly to the Ar case, we introduce the following notation for the generators corre¬ 
sponding to the roots: 




Cj C]^ , 





5 





2ej, 


(4.7) 
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Their forms are 


^(j-k) 

= a]ak - blbj, 



^{j,k) 

= Ojbk + albj, 

^(—j,—k) bjCik T bf^cij, 


^(hi) 

= 


(4.8) 


It is elementary to check the commutation relations, for example: 

[X(j _fc), X(A, _p] = [a]ak - blbj, alai - blbk] = a]ai - b\bj = -p, 

[X(j _fc), X(fc _j)] = a]aj - b]bj - a\ak + b\bk = Hj - Hk, etc. (4.9) 

The quadratic Casimir operator is 

r 

C 2 = Ntot{Ntot +‘^1^), Ntot =+ b^jbj), (4-10) 

i=i 

which gives M{M + 2r) in the present representation. It is easy to see that each number 
state belongs to some weight 

r 

|ni,..., ni,..., ^ /i = - nj)ej. (4.11) 

i=i 

In contradistinction with the A 2 r-i case this correspondence is not 1 to 1. Some weight 
spaces are degenerate. For example for M = 4 and r = 2, 


1 , 1 ; 1 , 1 ), | 2 , 0 ; 2 , 0 ), | 0 , 2 ; 0 , 2 ) 


belong to the null weight p = 0. 


As in the case of the binomial states (p.7|) we adopt as the ‘base’ state l-^o) the highest 
weight state 

(4)" 




'M 


| 0 >. 


(4.12) 


which guarantees “minimum uncertainty”. Together with all the positive root generators, it 
is also annihilated by the following generators: 


^{—j,—k)i hh/; 2 < j, fc < r, (4.13) 

which form a Cr-i subalgebra. Likewise the action of the Cartan subalgebra generator 
Hi does not change the highest weight state. Therefore the Cr multinomial states are 
parametrised by 

Sp{2r)/U{l) X Sp{2{r - 1)) = 


15 




which also indicates the connection to the A 2 r-i case. In fact the generators having non¬ 
trivial action on the highest weight state are 


2<j<r and I < j < r. (4.14) 

The generators in the hrst (second) group commute among themselves. In particular, 
which belongs to the lowest root, commutes with all the generators in the list 


(|4.14|) . The non-commuting pairs among the above generators are 

|V(-ij), V(_1 _,)] = -2A'(_1 _1), 2 <j<r, 


(4.16) 


and the resulting generator commutes with all the other generators in the list (|4.14|) , as 
shown above. 


In terms of 2r — 1 complex parameters 

0, 2 < j < r, e-.-, 1 < j < r, ^ = (6, .. .■ ■ ■,^-r) e 
the un-normalised coherent state is expressed as 


(4.16) 




|U/, - 

i=2 j=i 

with [C,D] = 2(X]j=2 -1) commuting with C and D. With the help of the B-C-H 

formula 

gC+D ^ gC-i[C,D]gD 

and the formal Taylor expansion theorem (|3.9|) we arrive at the following expression of the 
un-normalised Cr multinomial state 


M 


+ X! + X! I |0)’ 

i=2 i=i 


(4.18) 


in which the effects of non-commutativity cancel out exactly. Therefore the normalised Cr 
multinomial state is 


M 


|r/; M; Cr) = f ^ ri^a] + ^ r]_jb] | |0), 




(4.19) 


i=i 


in which 


“ 1 1 + XI l^-il 

i=2 j=i 


1 Vj=^jm, V-j=^-jVi, (4.20) 


satisfying the condition 

X(ihir+ih-iT) = 1- 

This has exactly the same form as the A 2 r-i multinomial state. 
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4.2 Probability Distribution 


Now we derive the probability distribution from the coherent state, which has exactly the 
same form as the multinomial state. So it predicts the multinomial distribution for the 
numbers ni,... ,nr with the corresponding probabilities |? 7 ip,... ,|?7_rp: 


|(ni, ...,nr;ni,... ,nr\r]-, M]Cr)\^ = 


M\ 




2nr.| |2ni 


Hi'. ■ ■ -Ur'-Ui'. ■ ■ -Ur'- 

As remarked above, the Cr states are labeled by the weight 


I ‘ 2 i 7 Xf 


(4.21) 


M (hi) • • •) hr) 


which takes positive, zero and negative integer values. Each weight space has one or many 
number states which are orthogonal to each other. Therefore the Cr multinomial distribution 
is obtained by summing the contributions from these number states: 


C^iri;M)= Y: 


M\ 




■Ihi 


|2ni 


Ihrl 1h-l 


2ni 


Ih-i 


|2nr 


(4.22) 


Let us interpret it in terms of “picking up balls from a pot”. The pot contains an inhnite 
number of balls of r-different colours. There are two types of balls for each colour, the 
^''positive” one and ''''negative'’^ one. Let the probabilities of picking one j-th colour ball be 
rfj for the '^positive" and for the "'negative”. We pick up total of M balls and ask the 
probability distribution for the “net” number of balls (or the “weight”) for each colour: 
lij = Uj — Uj, j = 1,... ,r. It is given by the Cr multinomial distribution. We see that the 
folding of the A 2 r-i Dynkin diagram leading to that of Cr is very suggestive of this situation. 


5 Dr Multinomial States 

Here we will derive probability distributions associated with the symmetric representations 
of Dr algebra. They have some new features not present in the multinomial distributions 
associated with A 2 r-i or Cr algebras. The Dynkin diagram of Dr algebra with the names of 
simple roots attached to the vertices is shown below. 



The corresponding simple roots are 

0(1 = 61 — 62, 0(2 = 62 — 63, . . . , O ' r -2 — (^ r -2 — Or-lj 0(^-1 = 6^-1 — Cr , 0 !r = 6^-1 + 6^. ( 5 . 1 ) 
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The positive roots are all of the same length: 


Cj-efc {j<k), Cj + ek- (5.2) 

The dimension of Dr algebra is 2r^ — r. The fundamental weights are 

Ai = Cl, A 2 = Cl + 02,..., (5-3) 

and we consider, as before, the irreducible representation with highest weight 


yU = M\i = Mei- 


(5.4) 


Let us denote this representation by and the corresponding vector space by . We 
know from Weyl’s dimension formula 


dim«) 


/M + 2r - 3\ 
V 2r - 3 ; 


X 


M + r — 1 
r — 1 


(5.5) 


Let us realise this representation in terms of 2r bosons 


Oi, . . . , Or, &1, . . . , 

and in its restricted Fock space denoted by F.^, 

|ni,... ,?7,r;hi,... ,nr), ni H- \-rir + ni-\ - \-nr = M. (5.6) 


We have 


dim(F“) 


/M + 2r - lA 
V 2r - 1 ) 


Comparing 


and ( b.7|) , we hnd 


/M + 2r - lA 

V M )■ 


dim{F^) = dim(V^) + dim(F 2 f-2) 

= dim(V^) + dim(V^-2) + ..., 


(5.7) 


(5.8) 


which means that the bosonic Fock space F^ contains several irreducible representations 
with different L’s. 

Let us introduce, as in the Cr case, the following notation for the generators corresponding 
to the roots: 


d^{j,k) T ^{—j,—k) ^k- (^"9) 
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Their forms are 


^{j,k) (^jbk ^k^ji ^(~jy~k) bjClk- 

It is elementary to check the commutation relations, for example they are 

[Xq-_ fc), X^k,i)] = Wjttk - blbj, alk - a}bk] = a]bi - ajbj = Xq- p 
[X(j-fc), X(_j _fc)] = a]aj - b]bj + a\ak - b\bk = Hj + Hk, etc. 

The quadratic Casimir operator is 

6-2 = {Ntot + 2 (r - 1 )) - 4XtX, N,,, = + b%), 

i=i 

in which K and are quadratic operators in the oscillators 

r M 

K* = J24>>'r 


(5.10) 


and: 


(5.11) 


(5.12) 


(5.13) 


i=i 


i=i 


They commute with all the above generators including those belonging to the Cartan sub¬ 
algebra: 

[K, X±(,-±,)] = [X, X,] = [Xt, X±(,-±fc)] = [Xt, X,] = 0. (5.14) 

In terms of X^ we can express the decomposition of the bosonic Fock space succinctly: 

= < © 1/^-2 © ... V^{V°), (5.15) 

in which the vector space is obtained from the highest weight state 

|M,0,...,0;0,...,0) = ^^|0), (5.16) 

by applying the negative weight generators successively. The j-th vector space in the right 
hand side is obtained from the highest weight state 

I' ©^-20-1) 

I iJ 


^(M - 2(j - 1))! 


(/V)^-‘| 0 ), 


by applying the negative weight generators successively. It is easy to see that X annihilates 


all the states in 


Kv = 0, Vn e V^, 

and we get C 2 = M{M + 2(r — 1)) in the highest weight representation ( b.4|) , (|5.16|) . 
easy to see that each number state belongs to some weight 

r 

|ni,..., ui,..., Ur) ^ /i = X!(X ~ 

i=i 


It i 


IS 


(5.18) 
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The highest weight state (|5.17| ) is annihilated by the following generators belonging to a 
Dr-i subalgebra 




Hj, 2<j,k<r, 


(5.19) 


as well as by all the positive root generators. The Cartan subalgebra generator Hi does not 
change the highest weight state. In other words, the generators having non-trivial action on 
the highest weight state are 


2<j<r. 


(5.20) 


If we denote the compact group corresponding to Dr by SO{2r), the Dr multinomial states 
are parametrised by 

SO{2r)/U{l) X ^0(2(r- 1)), 

having the dimension 

4(r — 1). 

In terms of 2 (r — 1) complex parameters 


0 ) 


we dehne a linear combination of the non-trivial generators (|5.20|) as 




(5.21) 


(5.22) 


It should be noted that all the generators in ( ^.22|) or (|5.2CI| ) commute among themselves, 
since the sum of the corresponding roots are not roots any more. Thus we arrive at the 
expression of the un-normalised coherent state: 


exp[T](al)^|0) = J] exp(0^(-ij)) H exp(e-iX(-i,-i))(al)^'^|0). 

i=2 j=2 


(5.23) 


By repeated use of the formal Taylor expansion theorem (p.9|) we obtain the following explicit 
form 


M 


J 


1 + ^- 3^3 ~ I |0)- 

j=2 j=2 j=2 


(5.24) 


This looks similar to the 242r-i and Cr multinomial states, except that the coefficient of b\ 
is not independent. The normalised Dr multinomial state is 


M 


\r];M;Dr) = 




(5.25) 
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in which 


Vi 


V-i 


(i + Eie.P + Eie-.P + lEe.e-.p) 

V i=2 j=2 j=2 j 

3=2 


Vj = ^jVu V-j = ^-jVi: 2 < j < r, 

(5.26) 


satisfying the condition 

E(i^ir+i^-iT) = 1- 

i=i 

Let us turn to the form of the probability distribution derived from the Dj. multinomial 
state, which has a form similar to that derived from the Ar multinomial state. Similar to 
the Cr case, Dj. multinomial state predicts the multinomial distribution to the number states 
with the probabilities and : 


\{ni ,... ,n^;ni,...,M; ^ = 


Ml 




Till • • • Tlrinil • • • Tlri 


2ni 


I Stt.t’ 


(5.27) 

By summing the contributions from all the number states belonging to a given weight /j, we 
obtain Dr multinomial distribution; 


D^(V;M)= -FT-(5.28) 

Thus the interpretation as “picking up coloured balls from a pot” is also valid. The marked 
difference is that among the probabilities ..., \Vr\‘^, • • • > only 2(r — 1) of 

them are independent. As is clear from ( b.26|) , one of the dependent probabilities, say |?7-ip, 
depends on the information of the other p±j’s including their phases (or more precisely 
^/s) , not |? 7 ±jps. We believe that this is a novel feature not encountered in any classical 
probability distributions. We may say that the Dr multinomial distribution has non-classical 
(or quantum) features. 


6 Br Multinomial States 

The Dynkin diagram of Br is obtained from that of Dr+i by folding the two tails. 
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Thus we expect that the Br multinomial states (distributions) have similarities with those 
of Dt with some added new features due to the folding. The simple roots of Br are 

Q-i = ei — 02, a2 = €-2 — 63, • ■ ■ , Q-r-l = Gr -1 ~ Cr; CHr — (6T) 

The positive roots are 

ej-Ck, {j<k), ej + ek, Cj. (6.2) 

There are 2r(r — 1) long roots and 2r short roots (±ej) and the dimension of Br algebra is 
2^2 _|_ sajne as xhe fundamental weights are 

Ai = Cl, A2 = Cl + 62, ... ( 6 - 3 ) 


As before we consider the irreducible representation with the highest weight 


H = MXi = Mci- (6.4) 

Let us denote this representation and the corresponding vector space by . Weyl’s 
dimension formula gives 


dim(y“) = 


+ 2r - 2\ 2M + 2r - 1 


2r - 2 


X 


2r - 1 


(6.5) 




This representation is realised in a restricted Fock space denoted by F’ 2 r+i 

, 77.J., 77-3, 


-^ 2 ^+ 1 ! |uo,ui,... ,?7,r;ui,... Uo + uiH-h + Ui H-h = M, 


( 6 . 6 ) 


which is generated by 2r + 1 bosonic oscillators 




61 ,. . . , 6 ^. 


As in the Dr case, by comparing the dimensions of the bosonic Fock space 


dim(F. 


M 5 
2r+l; 


+ 2r^ 
. 2r ; 


+ 2D 
. M . 


with the dimensions of (|6.5|), we find 


dim(F2f+i) = dim{V^) + dim{F^^^) 

= dim(V^) + dim(V^-2) + ..., 


(6.7) 


( 6 . 8 ) 


which means that the bosonic Fock space contains several irreducible representations 

with different highest weights {L = M, M — 2,... 


Pb 
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Similarly to the Ar case, the generators corresponding to various roots have the following 
forms: 

^{j,—k) ^k^ji 

^{j,k) ^{~jt~k) bj(^k ^k^ji 

^Ufi) cijCio ^i—j,—j) ttoOij bjCiQ, (6-9) 

in which, as in the Cr case, we use the notation: 




Cj , 










Cjj 

^(-i,o) -ej. 

(6.10) 


The commutation relations are easily verihed as in the previous cases. The quadratic Casimir 
operator is 

C 2 = Ntot{Ntot + 2r — 1) — AK^K, N^ot = (6-11) 

i=i 

in which K and are quadratic operators in the oscillators 

-I r 1 M 

^ = 9«o + II ^jbj, = ^(4)^ + £ (6-12) 

^ i=i ^ i=i 

As in the cases, K and commute with all the above generators including those 
belonging to the Cartan subalgebra. The decomposition of the restricted bosonic Fock space 
into the irreducible representation spaces goes in parallel with the Dr case: 


D 


M 

2r+l 




(6.13) 


in which the vector space is obtained from the highest weight state 

^(ol)"|0) = |0, M, 0 ,... i 0 ,..., 0), (6.14) 

by applying the negative root generators successively. The j-th vector space in the right 
hand side -g obtained from the highest weight state 


^(M - 2 U - 1))! 


(A-t)l->|0). 


(6.15) 


in a similar way. As in the Dr cases, K and annihilate all the states in . Thus 
the quadratic Casimir operator takes the value C 2 = M(M + 2r — 1) in the highest weight 
representation ( |6.4|) ,( |6.14| ). 
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One great difference between the Dr and Br cases is the correspondence between the 
number states and weights. In the Br case 

r 

|no, ni,..., Ur] ni,...,nr)^fx = XI (% “ ( 6 - 16 ) 

i=i 

Namely, hq, the number of ag quanta, has no effects on the weights. 

The Br coherent states can be constructed in a way similar to the Dr cases. The gener¬ 
ators having non-trivial action on the highest weight states are 

2<j<r, and X(_i,o), (6.17) 

which commute among themselves, since the sum of the corresponding roots are no longer 
roots. They constitute one half of the generators corresponding to the quotient space 

SO(2r + X SO{2r - 1), 


having the dimension 


In terms of 2r — 1 complex parameters 


2(2r - 1). 


^0, ^j, ^-j, 


(6.18) 


we dehne a linear combination of the non-trivial generators (|6.17|) as 


T - ^oX-1,0) + X + X 

j=2 j=2 


(6.19) 


Then the un-normalised coherent state is expressed as 


exp[T](al)^|0), 


( 6 . 20 ) 


which leads, after repeated use of the formal Taylor theorem (^.Oj), to 


+ X + X ~ (^ + X 1 |6)- 


M 


i=2 


t=2 


J=2 


( 6 . 21 ) 


Thus we obtain the normalised Br multinomial state 


M 


r];M;Br) = -j= r]ga\ + X + X |0)> 


'M\ 


i=i 




( 6 . 22 ) 
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in which 


i=2 


i=2 


C2 r 

I ^ 

2 


, Vo = ^oVi, 


i=2 


= |i + Elor + Eie-.r + l^ + Eoe-.i 

^2 r 

Vj = ijVu V-j = ^-jVi, 2 < j < r, T]_i = -(^ + 

^ j=2 


(6.23) 


satisfying the condition 


l^oT + E(I%I^ + = 1- 

i=i 


Let us turn to the probability distribution. The multinomial states give multinomial 
distribution to the number states with probabilities Ivol"^, IVjl"^ Iv-jl"^ ■ 


\{no,ni,... ,nr-,ni,... ,nr\r]; M; Br)\^ 
Ml 


(6.24) 


\V0\^^°\Vl\"^^---\Vr\"^^\V-l\"^^---\V-r 


I ‘2iTlr' 


riQinil • • • ririnil • • - rir 

By summing the contributions from all the number states belonging to a given weight /j, we 
obtain the Br multinomial distribution: 


BuimM) 


(6.25) 


= E 


Ml 






\v-^ 


I 2 ?t. 7 - 


Here let us recall that uq has no effects on the weights. Thus the interpretation as “picking up 
coloured balls from a pot” is also valid but with a slight modification. In the pot we have 2r+l 
types of balls, among them r different colours and each colour has “positive” and “negative” 
types. There are also “colourless” (or “dummy”) balls. They have probabilities Ivjl"^, Iv-jl"^ 
(j = 1,... ,r) and \vo\‘^- We pick up total of M balls and ask the probability distribution 
of the “net” number of coloured balls (or weights). It is given by the Br multinomial 
distribution. As in the Dr multinomial distribution, among the probabilities , 

1^-1 P, • • •, only 2r — 1 of them are independent. As is clear from (|6.23|) , one 

of the dependent probabilities, say \v-i\‘^, depends on the information of the other ri±/s 
including their phases. The existence of the “colourless” balls (or dummy elements) and the 
“quantum” nature of ? 7 _i are novel features of the Br multinomial distributions. 


7 Summary 

Starting from the fact established in our previous work [|T^ that the coherent states of 
the Heisenberg-Weyl, su{2), su{r + 1), sm( 1, 1) and su(r, 1) algebras in certain symmetric 
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(bosonic) representations give the well-known probability distributions, the Poisson, bino¬ 
mial, multinomial distributions with their “negative” counterparts, we have proceeded to 
the second stage in the study of “quantum probability”. By reversing the logic, we have 
obtained new probability distributions based on the coherent states of the classical algebras 
Cr and Dr in symmetric (bosonic) representations. These new probability distributions 
have similar features as the multinomial distributions related with Ar algebra. They also 
possess several new features reflecting their Lie algebraic and “quantum” backgrounds. As 
byproducts, simple proofs and interpretation of some addition theorems of Hermite polyno¬ 
mials are obtained based on the ‘coordinate’ representation of the (negative) multinomial 
states, the coherent states of su{r -|- 1) (sM(r, 1)) algebra in symmetric representations. 
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Appendix Addition Theorems II 


In this appendix we show a simple proof and interpretation of another type of addition theo¬ 
rems of Hermite polynomials. These theorems are non-compact counterparts of the theorems 
presented in section 3.3. They are obtained from the coordinate representation of the neg¬ 
ative binomial and negative multinomial states, i.e., the coherent states of the sm(1, 1) and 
su{r, 1) in symmetric representations. The theorem corresponding to the negative binomial 
states reads 


2_ (fcriZZ^iL 


(1 _ 

°° inl2Y 

~ X! -Ti- HnJrM-l{XQ)Hn{Xi)^ 


Xq - r]Xi 


n=0 


n\ 


(A.l) 


in which p is a complex parameter |? 7 | < 1. This addition theorem is known as generalised 
Mehler formula but is not found in the standard mathematics reference texts, except 

for the simplest case with M = 1 which is well-known as Mehler formula (see, for exam¬ 
ple, pl94 of |]^|). For a detailed characterisation of the negative binomial (multinomial) 
distributions in terms of Lie algebras, we refer to our previous work WA . 
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Let us begin with the negative binomial distribution (here 77 G R for simplicity): 

n = 0,l,..., (A.2) 


n 


which describes the probability distribution of the “waiting time” |^. Suppose we play 
Bernoulli’s trial of success and failure in which the probability of failure is 0 < < 1. The 

probability distribution for n, such that the (preset) M-th (M > 1, integer) success turns 
out at the M + n-th trial, is given by the above formula (|A.2D . We follow the examples 
of the previous sections and construct the “probability amplitude” of the negative binomial 
distribution. We choose the following restricted bosonic Fock space built by two bosonic 
oscillators: 


aj|0) = 0, j,/c = 0,l, 

ini 

\no;ni) = ° / , |0), no - ui = M - 1, n > 0. 


(A.3) 


Vnolnil 

Here no is the total number of trials except for the hnal one and ni is the number of failures 
(the hnal trial is always a success, by dehnition). Obviously this Fock space is inhnite 
dimensional. We look for a state \ri] M)~ such that 

For a special choice of the phases (cf. (|2.5| )) we arrive at a very simple result 

\V,M)~ = ^ |no;nl)(no;nl|?7;M)■ 

= (l-??^)^^|no;ni)p’’ 


no! 


ni!(M- 1)! 
— 0 (®o)^ ^ 


(A.4) 


This is called the negative binomial state |^, |^. This is exactly an sn(l, 1) coherent state 
as we will see presently. The sn(l, 1) algebra is realised in the above Fock space as 


A+ = a\a\, K_ = aoai, Kq = f-{No + Ni + 1), Nj = a]aj, 

2 

[iF+,iF_] = -2iFo, [iFo,iF±] = ±iF±. 


- 


The lowest weight state is \M — 1; 0): 


M, 


K.\M - 1; 0) = 0, Ko\M - 1; 0) = —\M - 1; 0), 


(A.5) 


(A.6) 
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which gives rise to the discrete irreducible representation with Bargman index M/2. Thus 
the un-normalised coherent state is (p G C) 


1;0) = 1;0), 

which has the same form as given in (|A.4|) . 


(A.7) 


Next we take the coordinate representation of the above negative binomial state: 




1 ; 0 ) 

and evaluate it in two different ways. The hrst is to simply expand the exponential and use 
the formula ( |3.33|) : 

^-\{xl+xl) ^ {r]/2y 


.,t„t 


(a:o;xi|e’'“o“i|M- 1;0) = 


E 


n\ 


■hfn+M—1 (^o)7fn(^l)) 


7rVy(M- 1)! ^0 

which corresponds to the right hand side of (|A.1|) . 

The second is to use the coordinate representation of the creation operators 
1, d . li„2„_i„2 „ d 


(A.8) 




to obtain 




(a:o;xi|e’'“o“i|M- 1;0) = 


- 1 ) 


M-l 


^^{xl+xl)^riDoDi/2 jjM-l^-{xl+xl) ^ 


(M — 1)! 

By applying the formal Taylor theorem (|3.9|) with respect to Xi by treating tiDq as a param¬ 
eter, we obtain 




1 ; 0 ) 

^ ^ ^ _ £)Ar-lg-(a;i+77Do/2)2g-x2 

7rVy(M- 1)! ° 






(A.9) 


7rVy(M- 1)! VI -h" 

which gives a scaled (l/\/l — V^) and shifted {—rjxi) Hermite polynomial (Hm-i) by Ro¬ 
drigues formula (^.341): 


r.h.s of (A.9) = 


7rVy(M- 1)! 

Here use is made of a simple formula 


2.-M. 1^2 

(1 — p ) 2 e2 og Hm-1 


Xq - TjXi 


(A.10) 































which can be proved, for example, by taking the Fourier transform. By comparing ( |A.9| ) 
and ( |A.10|) we arrive at the addition theorem of Hermite polynomials given above (|A.1|) . It 
should be remarked that the generalised Mehler formula ( |A.1| ) is also obtained from Mehler 
formula (M = 1) by differentiating with respect to Xq M — 1 times. 


Generalisation to the negative multinomial distribution 




2^M + ^1 +-h — l)!^2ni 


n!(M-l)! 

n = (no,ni,... ,n^), r/= (pi,...,G 
0 < = r]f-\ -h < 1, 




(All) 

(A12) 


is rather straightforward. We introduce a restricted Fock space generated by r +1 oscillators: 


%|0>=0. 


(A13) 


1^0) ■ ■ ■ 1 ^r) 


\nol 


4)' 


|0), no - (ni H- \-nr) = M-l. 


\/no\ni\ ■ ■ -rirl 

Then the “square root” of the negative multinomial distribution is 

)“ = G _ A/f — 1- 


\r] ; M) = (1 — T}^) 2 e 




(A. 14) 


which is an sn(r, 1) coherent state in an irreducible symmetric representation with the lowest 
weight state 

|M-1;0, ...,0). (A.15) 


The generators are 


K+j = ala], K_k = aoOfc, l<j,k<r, 
~ OjjO'k (j 7^ ^ 7^ 0); ly? ~ ajaj. 


(A.16) 


It is easy to see that they leave the combination 


A = No-{Ni + --- + Nr) 


and the above Fock space (A-131) invariant. Among the above generators the following r 
generators have non-trivial action on the lowest weight state (A-15|) 


K+j = ala], j = l,...,r. 

Thus in terms of r complex parameters rji,... , 1 ]^, satisfying the condition 

i^r = EiA-r<i, 


(A17) 


(A18) 
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we obtain an nn-normalised negative mnltinomial state 


1 ; 0 , ..., 0 ) = - 1 ; 0 , 


(A.19) 


which has the same form as (|A.14|) . By evalnating the coordinate representation of the above 
state ( A.19I) in two different ways, we obtain another form of addition theorem of Hermite 
polynomials: 


(1 - 772)-^/2e ° 


2 - Vr^rY 

" l-?)f- rif. 


Hm-1 


Xo - 


rjyXr 




Vr 


= E (A.20) 


n ,=0 


ni! 


Ur'. 


One can obtain this addition theorem by combining the addition theorems from the mnlti¬ 
nomial state (|3.32|) and that of the negative binomial state (|A.1|) , which reflects the fact that 
the negative mnltinomial state is also obtained by combining the negative binomial state 
and the mnltinomial state. 


Before closing Appendix, let ns mention another interesting form of addition theorems 
of Hermite polynomials which is obtained as a special case of ( |A.lj ). By setting xq = x and 
xi = 0, we obtain 

(1 _ = E (A,2i) 

V1 - T ^0 

Here nse is made of the relations 

H2n{0) = (-l)”(2n - 1)!! = (-l)n ■ 3 ■ ■ ■ (2n - 1), = 0. 

This form of addition theorems can also be obtained from another type of “coherent states” 
of su{l, 1) algebra. Let ns take the single boson Fock space (|2.6|) - (|2.8|) with the basis 
{|n), n = 0,l,...,} generated by a and ah The sm(1, 1) algebra is realised by 

/n = i(ay, A'_ = ia=, A'„ = lata + i (A.22) 

As before evalnate an nn-normalised “coherent state” 

e*-^+|M-1) = 1), |t| < 1, (A.23) 

in two different ways {t = —rf)- The above state is known as the ‘sqneezed nnmber state’ 
in qnantnm optics 0. for the ‘base state’ |M — 1) is not of lowest weight. 
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